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ABSTRACT 

It has been suggested that a resonance between a rotating bar and stars in the 
solar neighbourhood can produce the so called 'Hercules stream'. Recently, a second 
bar may have been identified in the Galactic centre, the so called 'long bar', which 
is longer and much flatter than the traditional Galactic bar, and has a similar mass. 
We looked at the dynamical effects of both bars, separately and together, on orbits 
of stars integrated backwards from local position and velocities, and a model of the 
Galactic potential which includes the bars directly. Both bars can produce Hercules 
like features, and allow us to measure the rotation rate of the bar(s). We measure a 
pattern speed, for both bars, of 1.87 ± 0.02 times the local circular frequency. This 
is on par with previous measurements for the Galactic bar, although we do adopt a 
slightly different Solar motion. Finally, we identify a new kinematic feature in local 
velocity space, caused by the long bar, which is tempting to identify with the high 
velocity 'Arcturus' stream. 
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1 INTRODUCTION 

The inner regions of the Milky Way contain a non- 
axisymmetric structure sometimes referred to as the tri- 
axial (or boxy) bulge and/or the Galactic bar. Primary 
evidence for this structure are the COBE/Diffyse In- 
frared Background Experiment (DIRBE) near infrared 
(NIR ) luminosity maps of the in ner Galaxy (|Binnev et al.l 
Il997l : iBissantz fe Gerhardl l2002h , but it has also been 
mapped via st ar counts using Two Micron All Sky Sur- 
vey (2MASS) (jLopez-Corredoira et al.ll2005h . 2MASS and 
Optical Gravitational Le nsing Experiment II (OGLE-II) 
(|Vanhollebeke et al.l l2QQ9h . Typical parameters for the tri- 
axial bulge/bar are axis ratios of 1.0:0.4:0.3, that it lies at 
an angle of between 10° and 40° to the line of sight to the 
Galactic centre (see e.g. table 1 of IVanhollebeke et al.[ l2009) 
and that its mass is of order 1 x 10 10 M , (IZhaol 1996; 
IWeiner fe Sellwoodlll999h . 

It has recently been proposed ([Benjamin et al.l l2005h 
that the inner Galaxy also contains a 'long bar', discovered 
as an over-density in star count data at 4.5 /xm taken in 
the Galactic Legacy Infrared Mid-Plane Survey Extraordi- 
naire ( GLIMPSE) surve y of th e Galactic centre made with 
Spitzer. iBeniamin et al.l (|2005h measure the half-length of 
this bar as 4.4 ± 0.5 kpc, and lying at a position angle of 
44 ± 10° to the line of sight to the Galactic centre. Only 
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the closer end of the bar was detected, as the far end would 
lie at apparent magnitudes fainter than the completeness 
limit of the survey. The long bar has been confirmed by 
lLopez-Corredoira et al.l ([2007T L who detect it via red clump 
stars identified from the 2MASS survey. They find a half- 
leng th of 3.9 kpc and a p osition angle of 43 ±7°, very similar 
to IBeniamin et al.l (|2005h . The long bar is much longer than 
any previously proposed bar/bulge in the Galactic centre, 
for which the long ax is is thought to lie in the range 2.4 
(e.g iDwek et alJll995h to 3.5 kpc (e.g. IBissantz fe Gerhardl 
l2002n . Most interestingly, the long bar seems to lie at a dif- 
ferent angle to the traditional bar/bulge (hereafter 'Galactic 
bar'), fo r which estimates typicall y lie in the range 25 =b 10 
degrees (jVanhollebeke et al.ll2009h . 

Most recently, ICabrera-Lavers et al.l (|2007h have used 
wider field 2MASS data to isolate red clump stars over a 
wide range of Galactic longitude, showing that the long bar, 
the triaxial bulge and the disc can all be seen as separate 
structures (at least on our side of the Galaxy) with the im- 
portant result that the long bar is much flatter than the 
bulge. The long bar is found to have a vertical scale- height 
of around 100 pc whereas the bulge has a vertical scale- 
height of 500 pc. Furthermore, the red clump stars used to 
trace the long bar can be used to estimate its total mass 
(assuming that the mass to light ratio of red clump stars is 
the same in the long bar as it is in the bulge) as 6 x 10 9 
M©, of th e same order of the mass of the Gal actic bar 1 x 
10 10 M , (|Zhaolll996l : IWeiner fe Sellwoodlll999h . 
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Figure 1. Local velocity space, data includes all the stars 
with U and V velocit ies f rom the Geneva- Cope nhagen survey 
(iHolmberg et al.ll2007h and ISchuster et al.l (2006). The Hercules 
stream is the prominent feature area around (—25,-30) km s — 1 . 
The velocities have b een corrected to the solar velocity from 
ISchonrich et alJ (|2010h 

iDehnenl (2000) has shown that a bar in the inner galaxy 
has interesting dynamical consequences for the orbits of 
stars even in the solar neighbourhood, 8 kpc from the centre, 
due to a resonance of stellar orbits with the Outer Lindblad 
Resonance (OLR) of the bar. The manifestation of t his res- 
onan ce in our local space is the £ Herculis stream (|Eggenl 

Il971bh , or the 'Hercules stream' as it is usually called . 

In this paper, we use a similar method to IDehnenl 

(2000), to reexamine the dynamical effects on local disc 
stars of these non-axisymmetric features in the Galac- 
tic centre. We have e xamined both the long bar (see 
Lopez- Corredoira et al. | | 2007f ) and the Galactic bar (see 
Bissantz &; Gerhardll2002l V separately, and as a dual-bar sys- 
tem. We adopt R© = 8 kp c throu g hout the paper. This fa- 
cilitates comparison with IDehnenl $2000), who adopted the 
same value. 

We have made a probability-statistical study of the ef- 
fects of the bar(s) on local velocity space, by the means of 
potential theory. We will present the statistical evidence pro- 
duced by various variations of the long bar-system, and a the 
Galactic bar- system, as well as, some interesting combined 
models. 



2 THE HERCULES STREAM 

In Fig. [T] we show local velocity space in (£7, V), where U is 
the velocity towards the Galactic centre and V is the veloc- 
ity in the direction of Galactic rotation. The figure has been 
prepare d using 12939 s tars from the Geneva Copenhagen 
Survey (|Holmberg et al.f2007h of abundances, ages and kine- 
matics of nearby stars, and has been supplem ented by 1535 
stars f rom the the high velocity catalogue of ISchuster et all 
(2006). The bin size is 2 km s _1 , and the contours are 



in 5 percent steps relative to the maximum density. The 
veloc ities have been c orrec ted for the solar motion taken 
from ISchonrich et alJ (|2010h of (Uq,Vq) = (11.1,12.2) km 
s _1 . The Hercules stream is the feature at approximately 
(U,V) = (- 25,-30) km s" 1 . This is the feature which 
IDehnenl (2000) identifies as having a dynamical origin as 
a resonance with the bar. Our simulations in this paper of a 
Galactic bar, long bar (and both) are primarily constrained 
by this feature. 

The bulk of the stars in Fig. [1] belong to the compli- 
cated feature above the Hercules stream. The classically 
identified features in this region are the Pleiades, Sirius and 
Hyades streams, which occur at around (0,-10), (20,20), 
and (—25,-5) km s _1 respectively. Further streams than 
those mentioned have been identified by a number of groups, 
but these mainly concern higher velocity stars and cannot 
be pointed to easily in Fig[T] One of the very interesting fea- 
tures at higher velocities is the Arcturus stream at V ~ 100 
km s~* (for a wide range of U velocit ies) identified b y Eggenl 
(|l971ah and discussed recently by Williams et al. (2009), 
who have identified it in the RAdial Velocity Experiment 
(RAVE) survey. The Arcturus stream may have a dynami- 
cal origin, as we show using the simulations in this paper. 



3 MODELLING THE DISC AND BAR 

We model the effects of a bar on local disc stars by setting 
up a potential describing the disc, bulge, dark halo and bar, 
and integrating a library of orbits pa ssing through the Solar 
neighbourhood, similarly to IDehnenl (2000). 

We beg i n with the axisymmetric Galactic potential of 
iFlvnn et al.l (p.996). This model contains a disc, bulge and 
dark halo. The scale-length of disc matter (Hr) in the model 
is 4.4 kpc, which these days appears rather long, as recent 
studies of the disc scale-length give results closer to 3 kpc. 
Analysis of NIR data in J and K with DEep Near Infrared 
Surv ey of the Southern Sky (DENIS) give h R = 2.3 ± 0.1 
kpc (jRuphv et al.| [l996). while K giants s een in 2MASS give 
Hr = 3.3 ±0.1 kpc (|L6pez-Corredoira et"aD 120021 ). These 
surveys tend to probe the scale-length of luminosity in the 
disc, rather than the mass, which is traced rather better via 
main sequence dwarfs. Hubble Space Telescope (HST) star 
counts of M d warfs yield a disc scale- length of Hr — 3.0 ±0.4 
(Gould et al I ll997h . F and K dwarfs, traced in huge num- 
bers in the 6500 de g 2 SPSS survey, y ield a disc scale-length 
of h R = 2.6 ± 0.5 (Juric et al. 2008). For a more complete 
review of the disc scale- length, see iGardner et al.l (|2QQ9h . 
We chose a disc scale-length of 3 kpc to reflect the com- 
bined observational constraints. This was straightforward 
and only involved modifying the Miyamoto discs from which 
the exponential- like disc is built up. 

The surface density (£©) of the new disc at the 
Sun's position is 50 M© pc -2 and the local density p© 
i s 0.10 M© pc~ 3 , consist ent with observational constraints 
(|Holmberg &; Flvrmll2000h . The surface density of the disc as 
a function of Galactocentric radius is shown in Fig. [2] The 
new disc model has a 3 kpc scale-length, with a truncation 
at about 18 kpc. 

The full parameters of the new disc model are shown 
in Table [T] along with the other parameters of the model 
(which are unchanged but are reproduced for clarity). 
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Table 2. Default simulation parameters. 



10000 15000 

R/pc 



20000 



25000 



Parameter 



Value 



Unit 



Galactic Bar 

Angle 
Dimensions 

Mass 

Pattern Speed 

Local standard of rest 

OLR 



25 

3.5:1.4:1.0 

10 

55.9 

239 

1.87 



kpc 
10 9 M 
km s — 1 kpc -1 
km s — 1 



Long Bar 

Angle 

Dimensions 

Mass 

Pattern Speed 

Local standard of rest 

OLR 



43 

3.9:0.6:0.1 
6 

54.9 
235 
1.87 



kpc 
10 9 M 
km s — 1 kpc -1 
km s — 1 



Figure 2. The surface density of the disc component (solid line) 
as a function of Galactocentric radius. The dashed line corre- 
sponds to an exponential density falloff, 3 kpc, which is a good 
fit to the model over a wide range of radii. Note that the density 
truncates strongly at 18 kpc. 



Table 1. Parameters of the model. 



Parameter 


Value 


Unit 


Vh 


220 


km s _1 


ro 


8.5 


kpc 


M Cl 


3 


10 9 M 


r Cl 


2.7 


kpc 


Mc 2 


16 


10 9 M 


r c 2 


0.42 


kpc 


M dl 


77.04 


10 9 M 


adi 


5.81 


kpc 


M d2 


-68.48 


10 9 M 


a d 2 


17.43 


kpc 


M d3 


26.75 


10 9 M 




34.84 


kpc 


b 


0.3 


kpc 



The equations for the axisymmetric potential are: 



$ = $H + $C + ®D + $bar 



d> H = Iv^ln(r 2 +r 2 ) 



GM Cl 



GMc 2 



2 

C 2 



-GM Dn 

^j{R 2 + {a dn + ^{z 2 + b 2 )) 2 ) 



2 , 2 , 2 

x +y + z , 



(i) 

(2) 
(3) 

(4) 



and R = x 2 + y 2 . 



where = E„ =1 $D n , r = 

We have used a different m etho d to modelling t he ba r 
as was adopted in IDehnenl ((2000) and lMinchev et al . (2007). 
Instead of modelling the bar as a local quadrupole perturba- 
tion to the potential, we model the bar directly in the inner 
galaxy. This means that orbits of stars entering into the in- 



ner disc regions and even passing through the bar are much 
more accurately modelled, even if it is computationally far 
more expensive. 

T he bar is m o delled as a Ferrers n — 2 potential, as laid 
out bv lPfennigerl (|l984h . and as such, is a triaxial ellipsoid. 
The equations for t he bar are too c omplex to present here 
(for full details see |Pfennigedfl984h . The bar rotates as a 
rigid body. 

A bar is specified by its three axis lengths, its mass, its 
angular (or pattern) speed and its angle to the line connect- 
ing the Galactic centre and the Sun. 

We hav e two standard bars , with parameters 
taken from Lopez-Corredoira et al.l (|2007l ) and from 
iBissantz fe Gerhardl (|2002h . The dimensions of the 'long bar' 
are (3.9:0.6:0.1) kpc, with a mass of 6 10 9 M and an an- 
gle of 43°. The 'Galactic bar' has dimensions of (3.5:1.4:1.0) 
kpc, a mass of 1 10 M a nd an angle of 25°. These masses 
correspond to a iDehnenl (|2000h a value of 0.0037, for the 
long bar case, and 0.0040, for the Galactic bar case. Note 
that in the simulations we do not vary the dimensions of the 
respective bars, only the mass and position angle. 

An impo rtant paramete r is the pattern speed of the 
bar. Following IDehnenl (|2Q0Qh and lMinchev et al.l (|2007 ). we 
parametrise the bar's pattern speed by the OLR, i.e. by 
the ratio of the period of a star with velocity of the local 
standard of rest (LSR) in a given Galactic model to the 
period of rotation of the bar. 

The definition of the LSR is straightforward in axisym- 
metric models of the Galaxy, but needs to be redefined 
slightly once we have added a bar to the model. In the pres- 
ence of a bar, a star can no longer follow a circular orbit at 
the solar radius, but rather oscillates around a mean radius. 
For typical bar parameters, the oscillation is up to 10 km 
s _1 , compared to a mean velocity of 220 km s _1 in typical 
models. For each model run in our simulations, we deter- 
mine, experimentally, at what velocity a star, which starts 
at the Solar position, would have a mean orbital radius of 8 
kpc. Integration is carried out for 10 Gyr, so that the effect 
of the effect of the angle of the bar is averaged out. This is 
denned as the LSR for that model. 
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4 SIMULATIONS AND METHODS 

The basic simulation approach is to model the appe arance of 
local velocity space using the method described by iDehnenl 
(2000), a 'backward restricted iV-body method'. The idea is 
to trace a library of orbits starting from the solar neighbour- 
hood backwards over some fixed time (typically 1 Gyr) , and 
use the known density distribution and velocity distribution 
of the disc to reconstruct the velocity distribution of stars 
in the solar neighbourhood. 

We introduce two refinements, firstly, a more realistic 
representation of the bar potential and secondly, we use ob- 
servational constraints on the velocity dispersion and asym- 
metric drift of stars in the Galacti c disc over a wide range o f 
Galactocentric radius obtained bv lLewis &; Freeman! (|l989h . 
In practise, this latter refineme n t is n ot very different to the 
approach adopted by IDehnenl (|2000h , (simple exponential 
functions of stellar density and radial velocity dispersion for 
the old stellar disc) and indeed we find that we reproduce 
his results quite well when we adopt a bar simila r to his. 
Our simulations are in two dimensions (2D), as in IDehnenl 
((2003). 

The initial conditions are set by a velocity- grid, of 100 
by 100 elements in the UV- velocity plane, covering U — 
-100 to 100 km s" 1 , and V= -150 km s" 1 ) to V = 50 km 
s _1 ), in steps of 2 km s _1 . A tracer with the UV- velocity 
selected from the grid is then set into motion, following the 
tracer backwards in time over a period of 1 Gyr, tracking 
its position every 1 Myr. The initial position is selected as 
the solar position (8 kpc), and both the z-coordinate and 
the W- velocity are set to zero. The initial velocity in the 
rotational direction V is, of course, relative to the LSR for 
particular model being tested. 

Statistically, we use a density and velocity-based 
system, where a weight value corresponds to a 8 kpc- 
centric exponential distribution, and a Gaussian velocity- 
distribution, with sta ndard deviations as described in 
Lewis & Freema The V-velocity is corrected back 

to a 'true' V- velo city using the LSR of th e position and the 
asymmetric drift ([Lewis fe Freemad[l989h . 



5 RESULTS 

We examine in this section the effects of both the 'Galactic 
bar' and the 'long bar' on local velocity space. We vary the 
bar mass, value of the OLR (i.e. the pattern speed of the 
bar), and the angle of the bar. Finally, we look at the effects 
of including both bars at the same time. 

5.1 Effect of the mass of the Galactic bar 

We begin by running the standard model of the 'Galactic 
bar' but varying its mass over a range of 1 to 10 x 10 9 
M , showing the results in Fig. [3] The OLR for this case is 
1.87. The variations of mass within one orde r of m agnitude 
have prominent effects. As found bv IDehnenl (|2000h we find 
that a Hercules- like feature arises in local velocity space in 
the presence of a bar. We will henceforth refer to this as a 
secondary feature in velocity space, with the primary fea- 
ture (above it) containing most of the stars. The effect of 



increasing the mass of the bar is similar to that found by 
IDehnenl (|2000h , with the strength of the secondary feature 
increasing (and the depth of the resonance gap between the 
secondary feature and primary feature also growing). The 
secondary component becomes more ellipsoidal, as mass is 
increased, and, as expected, becomes more prominent. Also, 
the position of the secondary feature moves to more negative 
V velocities as the mass is increased. While it is clear that 
something like the Hercules feature can be generated by this 
bar, the detailed match is not yet perfect. In Fig.[T] the Her- 
cules feature is found at (£7, V) — (—25, —30) km s _1 . In the 
simulations, the generated feature approaches this location 
for the higher mass models, but at the expense of generating 
a sharp 'tail' of stars at positive U values spreading outward 
from the pr i mary featur e. This tail of s tars i s also seen in 
the IDehnenl {2000) and iMinchev et~ai1 (|20Q7h simulations, 
and increases with the mass of the adopted bar. Such a tail, 
if present at all in the observational data (Fig. [TJ is quite 
weak. Apart from this tail, the generated feature is fairly 
consistent with various estimates of bar mass where they 
are consistentl y around and o ver 10 10 M see e.g. Zhao 
Il996l : IWeiner fe Selrwoodlll999h . 

5.2 Long bar 

We next look at the 'long bar' case. We adopt the standard 
model for its dimensions and initial position angle and vary 
only the mass. These models all give a Hercules-like sec- 
ondary feature, similar to the Galactic bar case. The effect 
of varying its mass on local velocity space is shown in Fig. 2] 
One can see that mass takes a significant role in the division 
of the primary and secondary feature, where the secondary 
is stronger with increasing mass. Again, these plots superfi- 
cially resemble the real local velocity space (Fig. [TJ but not 
in detail. The position of the secondary feature tends to be 
at too high V velocities and the tail of stars at high U values 
also appears for the higher mass models. 

One very interesting 'tertiary' feature appears with the 
long bar, located at V — —80 to —100 km s _1 . The tertiary 
feature becomes more prominent as mass increases. This fea- 
ture is not seen in the Galactic bar simulations (Fig. [3]). We 
will return to this feature in section 15.41 



5.3 Pattern speed of the bar 

As shown bv IDehnenl (2000) and IMinchev et al ] (|20Q7h . the 
position of the Hercules feature is quite dependent on the 
adopted pattern speed of the bar, parametrised by the OLR. 
The feature can therefore be us ed to measure the pattern 
speed of the bar. IDehnenl |2000) finds that the feature can 
be fit by an OLR of 1 . 85 ± 0.15, using similar simulations to 
ours. IMinchev et al.1 (|2007h find a similar and more tightly 
constrained value of the OLR of 1.87 =b 0.04, using similar 
simulations as here but also from constraints implied by the 
locally measured Oort constant C. 

Our standard model for both the Galactic and long bars 
assumes an OLR of 1.87. For the Galactic bar case, the effect 
of shifting from an OLR value of 1.87 to 1.90 is shown in Fig 
\E\ The secondary feature shifts by about 5 km s _1 to more 
negative V velocities. A more detailed view of this is shown 
for the long bar case in Fig. [6] where we have varied the 
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Figure 3. Effect on local velocity space of varying the mass of 
the Galactic bar, from 10 9 to 10 10 M@ in increments of 10 9 Mq. 
For other simulation parameters see Table [2] 

value of the OLR ranging from 1.83 to 1.95. The secondary 
feature shifts quite markedly away from the primary feature 
over this range of assumed OLR values. We derive a best 
fitting OLR of 1.87 =b 0.02, assuming that the position of the 
secondary feature is at V — —30 =b 5 km s _1 . 

Our determin ation of the OLR value for the b a r is in 
agreement with iDehnenl |2000) and iMinchev etHI (|2007T ) 
within the error bars. We note that our value about the 




Figure 4. Variations of velocity space with different values of for 
long bar mass. Masses are 4, 5, 6 (middle), 7, and 8 1O 9 M0. For 
other simulation parameters see Table [2] 




Figure 5. Variation of velocity space with a Galactic bar model, 
OLR values of 1.87 and 1.90. 

same as theirs, even t hough we adopt the solar motion of 
ISchonrich et al l (|2010h . which shifts this feature by 7 km s 
to more positive V velocities in the observational plane (the 
Hercules feature is centred o n about V = — 35 km s _1 with 
the solar motion adopted bv lDehnenl (|2000h . rather than at 
V = —30 km s _1 as in our Fig.[T] Had we adopted the older 
solar motion, we would obtain an OLR of 1.90 =b 0.03. 

5.4 Bar angle 

IDehnenl (|2000h used the Hercules stream to probe the angle 
of the bar, deriving an angle of 25°. This is to be compared 
to the observational evidence for the (Galactic bar) angle, 
which lies in the range 20° to 40° . We probe the effect of bar 
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Figure 6. Effect on local velocity space with different values of 
for the OLR for the long bar case. Values are 1.83, 1.85, 1.87, 1.90, 
1.93, and 1.95, from left to right, top to bottom. There is a clear 
shift in the position of the secondary feature to lower V velocity 
as the OLR increases, allowing a determination to be made of the 
pattern speed of the bar. 

angle as well, for both the Galactic bar and long bar cases. 
As found by iDehnenl ([2000h , the effect is to change of the 
position of the secondary feature to towards more negative 
values of U and more positive values of V with increasing 
bar angle, and also to change the strength of the feature, as 
seen in Figs [3 and Fig. [8] However, these changes are rather 
subtle and we consider them too weak to constrain the bar 
angle. 

The tertiary feature is clearly dependant on bar angle, 
both for the long bar case and the Galactic bar case. The 
feature appears in both cases when the initial bar angle is 
greater than 40° . The orbits of stars in this feature are highly 
eccentric and take them to within 2 kpc of the Galactic cen- 
tre, ideal for direct interaction with a bar. This accounts for 
why bar angle is important to the generation of the tertiary 
feature. 

There is a stream in the Solar neighbourhood at about 
this ve locity, V ~ —100 km s _1 , the Arcturus stream ([Eggenl 
Il971ah . T he stream has re c ently b een recovered in the RA VE 
survey bv IWilliams et all (|2009h . IWilliams et~aH (|2009h ar- 
gue against an accretion origin for this stream, because its 
stars have a wide range of metallicity, so that the accreted 
system would have to have an unreasonably high mass, per- 
haps as large as the Large Magellanic Cloud (LMC). It is 
interesting that our simulations show that a dynamical ori- 



Figure 7. Variation of Galactic bar angle, 15°, 25°, 35°, and 45° 




Figure 8. Variation of long bar angle, 23°, 33°, 43°, and 53°. 

gin for the stream is plausible. Streams with a dynamical 
origin are expected to have a wide range of metallicity, as 
wide as the stars of its parent population. If so, the more 
plausible bar for originating this stream is the long bar, be- 
cause its measured bar angle is 43°. The Galactic bar could 
also produce this stream, but it would have to also be at 
an angle of about 45° , whereas most studies find that it lies 
much close r to the line of sight to the Galactic centre, at 
about 20° (jVanhollebeke et al]l2009h . 

5.5 Integration time 

We have chosen and integration time of 1 Gyr, but this 
choice is somewhat arbitrary. Fig. [9] shows the effects of 
adopting shorter and longer integration times of 500 Myr 
and 2 Gyr for the long bar. 
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Figure 9. Variation in integration time, 500 Myr, and 2 Gyr, 
using the long bar. 

After 500 Myr, the secondary feature has started form- 
ing but is not yet prominent. This represents about 4 bar 
rotations, and it is not surprising that it takes some time 
for the resonance with the bar to develop the trough be- 
tween the primary and secondary features. After 2 Gyr we 
can see the stronger splitting of the primary and secondary 
feature, with the resonance- area clearing, as well as the pri- 
mary feature tailing out to more positive U. Simulations as 
long as this assume of course that the bar has existed for 2 
Gyr with the same parameters, in particular that the rota- 
tion rate and mass have not changed during this time. Our 
choice of 1 Gyr is meant to allow structure enough time to 
develop, while not taxing the credibility of bar prope rties 
being very long lived in galaxies like the Milky Way. Fux 
(2001) has pointed out a potential difficulty with the back- 
wards integration technique used here. He finds that as one 
integrates further backwards in time, increasingly spurious 
structures develop in the UV -plane. We find no such be- 
haviour over approximately 16 periods of the bar (Fig. [9}. 
We are confident in our 1 Gyr timescale, for interpreting the 
simulations. 




Figure 10. Dual bar models, phase-locked 1.87 and 1.90 (top), 
unlocked, long at 1.87, galactic at 1.5 (middle left), long at 1.5 
galactic at 1.87 (middle right), half-mass locked at 1.87 (bottom 
left), half-mass unlocked long at 1.87, galactic at 1.5 (bottom 
right). 



5.6 Dual bar models 

Our simulations have shown that both the Galactic bar, and 
the newly discovered 'long bar' can generate features in ve- 
locity space which are superficially similar to those seen in 
the Solar neighbourhood, in particular the Hercules stream. 
We now investigate the effect of including both bars. 

We show a selection of dual bar models in Fig. 1101 We 
begin with both bars having the same OLR - i.e. the two bars 
have the same pattern speed and rotate together. The up- 
per panels of Fig. HOIshow the local effect on phase space for 
OLR values of 1.87 and 1.90. The overwhelming impression 
here is that the main feature in phase space is driven to large 
values of U, inconsistent with observations. The main cause 
for this is simply the large amount of mass which now resides 
in the Galactic centre, in these cases 1.6 10 10 M . Similar 
behaviour is seen even for cases with a single bar, when the 
mass of the bar exceeds about 10 10 M©. More than this 
much mass and the bar(s) starts to generate an unaccept- 
able amount of activity on the Solar neighbourhood (£/, V) 
distribution, typically pushing stars to too large values of U. 

Two decoupled scenarios were tested, with one of the 
bars at an OLR of 1.50 and the other at 1.87. These are 
shown in the middle panels of Fig. [10] On the right, the 
Galactic bar has an OLR of 1.50, while the long bar has 
an OLR of 1.50 on the left. Again, the total mass in the 



two bars is generating significant structure in local veloc- 
ity space, with the same problem of too many stars being 
shifted into a tail of high U velocities. In addition, the ter- 
tiary feature at V ~ — 80 becomes very prominent, as the 
chances of an interaction with either bar is now higher. We 
tested this by running single-bar models with an OLR of 
1.5, which results in no secondary or tertiary feature for 
the Galactic bar case, and in the long bar case it causes no 
secondary and a very weak tertiary, such as seen in other 
single-bar simulations of the long bar. Both the upper and 
middle panels show more structure than really exists in the 
Solar neighbourhood, primarily because too much mass is 
assigned to the bar structures in the inner Galaxy. To alle- 
viate this problem, we tested a scenario where the masses 
of both bars are simply halved. This is shown in the lower 
panels of Fig. 1101 On the left, the bars both have an OLR 
of 1.87, while on the right, the bars are unlocked, with the 
Galactic bar at an OLR of 1.5 and the long bar with an 
OLR of 1.87. Reducing the mass brings the simulations into 
much better agreement with the observations, although the 
secondary feature is at a too high V velocity. This can be 
remedied by shifting to a higher value of 1.90 for the OLR. 
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6 DISCUSSION & CONCLUSIONS 

We have performed simulations of the effects of bars in the 
inner Galaxy on the distribution of stars in veloci ty space 
near t he S un. Our simulations ar e similar to those of lDehnenl 
(|2000h and lMinchev et all (|200 /t ). in which a large library of 
orbits of stars passing through the Solar neighbourhood are 
performed in a model of the Galactic potential including a 
bar. These earlier studies showed that the Hercules stellar 
stream in the solar neighbourhood could be the result of a 
resonance of stellar orbits with a fast moving bar (i.e. with a 
co-rotation with t he Galactic dis c not much beyond the end 
of the bar). Both [Pehnenl ([2000) and iMinchev et al.1 (|200/t ) 
find that the Hercules stream can be used to constrain the 
rotation rate and angle of the bar. 

Our motivation to reexamine this is the recent discovery 
that the Galaxy may contain a second bar, in additio n to 
the traditional Galactic bar, b y lBeniamin et al.l (2005) and 
lLopez-Corredoira et al.l (|2007h . These authors have found a 
'long bar' in the Galactic centre. It is quite different to the 
Galactic bar, being highly flattened (semi-major axes of 3.9, 
0.6 and 0.1 kpc as compared to 3.5, 1.4 and 1.0 kpc). Fur- 
thermore, it is aligned at an angle of about 43 degrees to the 
line of sight to the Galactic centre, as opposed to about 20 
degrees for the Ga lactic bar. For the Galactic bar we ad opt a 
mass of 10 10 M (|Zhaolll996l : [Weiner fe Sellwoodll 19991) an d 
for the long bar 6 x 10 9 M (jLopez-Corredoira et al.ll2007l ). 
The long bar has a mass quite comparable to the traditional 
bar. 

We have simulated the effects of the long bar, as well 
as the Galactic bar, and both bars, on the velocities of stars 
in the Solar neighbourhood. Rather than simulating bars as 
a quadrupole perturbation in the local Galactic potential, 
we simulate the bars with a Ferrers potential (which gener- 
ates a triaxial ellipsoidal density distribution) . This is com- 
putationally more expensive, but allows for more accurate 
modelling of a bar as a whole, especially for stars passing 
through, or trapped in, a bar. The simulations are performed 
in 2-D, although the model itself allows 3-D simulations. 

We confirm the basic picture that the Hercules stream 
can be generated by a resonance between local stars and 
the bar. Both the long bar and the Galactic bar produce 
Hercules- like features in the Solar neighbourhood. The po- 
sition of the Hercules stream in the simulations is found to 
depend on the mass of the bar, the rotation rate of the bar, 
and rather weakly, on the angle of the bar. The geometry 
of the bars is not very important either, since both bars 
generate Hercules-like streams. 

We measure the rotation rate of the bars from the po- 
sition of the Hercules stream, finding that both bars can 
produce acceptable fits to of the (U, V) velocities of nearby 
stars. We measure the rotation rate of both bars via the 
Outer Lindblad Resonance value they produce (i.e. by the 
ratio of the period of a star with velocity of the Local Stan- 
dard of Rest in a given Galactic model to the period of rota- 
tion of the bar). We measure values of 1.87±0.02 for the OLR 
of both the Galactic bar and the long bar. This value gives 
a pattern speed of 55.9 and 54.9 km s _1 kpc -1 , respectively, 
for the Galactic bar and long bar. This puts the correspond- 
ing co-rotation, for the long bar, as defined by its potential 
model, at the approximate tip of the bar, 3.9 kpc. For the 
Galactic bar, our model produces no co-rotation radius. It 



should be noted however that this value is sensitive to the 
assumed S olar motion. We a ssume the Solar motion recently 
derived bv lSchonrich et al.1 (l2010h. Had we adopted the So- 
lar motion used bv iDehnenl (|2000l ). we would have derived 
a value of 1.90 =b 0.03. Bot h of these deri ved values are in 
reasonable agreem ent with jDehnen (l2000l Vs determination 
of 1.85 ± 0.15 and lMinchev et al.l (|2007l Vs determination of 
1.87 ±0.04. 

After looking at the effects of each bar alone, we tried 
putting both bars into the modelling. The effects on local 
velocity space are quite dramatic, as we have now added 
quite a lot of mass in an non-axisymmetric component to 
the Galactic centre. In particular, these simulations produce 
a striking tail of stars in the Solar neighbourhood at high 
U velocities, inconsistent with observations. The simple ex- 
pedient of halving the mass of both bars produces much 
better fits to the data. On the basis of these experiments 
then, if there are two bars in the Galactic centre, we won- 
der whether the mass estimates of each bar include material 
from the other bar, and are perhaps overestimates. This of- 
fers a challenge to observational studies of the bar masses. 

Finally, we have found that an additional feature in lo- 
cal velocity space can arise under some conditions, especially 
for the long bar case. A weak stream of stars can form at 



about V 



-100 and U < km s . We believe this is 



due to a direct interaction of such stars with the bar, rather 
than a resonance with the bar. It is tempting to associate 
this feat ure with the known Arcturus stream at —100 
km s -1 ( Willi ams et al.l [20091 ) in the Solar neighbourhood. 
Thus, a dynamical origin for the stream is possible, whereas 
people have mainly discussed accretion origins (of a satellite 
galaxy) to date. 

We point out that our method has shortcomings. 
Firstly, the simulations are in 2-D, and while this is almost 
certainly acceptable for typical disc stars with near circular 
orbits, full 3-D simulations may be needed for features like 
that which we associate with the Arcturus stream. In ad- 
dition, there are no spiral arms or other non-axisymmetric 
components, besides the bars, in the modelling, and these 
may be responsible for additional complexity in local veloc- 
ity space which we cannot model. We are not affected by the 
problems of long-term backwards integration as mentioned 
m lFul (|200lh . Compared to earlier studies, we try to more 
accurately represent the Galaxy's observed properties, such 
as local densities and bar geometry. The model also lets us 
look at direct interaction, as well as resonances, this is es- 
pecially important for high velocity stars (e.g. those in the 
Arcturus stream), where they have high enough eccentrici- 
ties to enter the inner parts of the Galaxy. The model also 
makes it possible, in the future, to fully study the 3D impact 
of the bar(s). 
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